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ABSTRACT 

We describe a pseudo-3D photoionization code, NEBU_3D and its associated visualiza¬ 
tion tool, VISJS1EB3D, which are able to easily and rapidly treat a wide variety of nebular 
geometries, by combining models obtained with a ID photoionization code. The only re¬ 
quirement for the code to work is that the ionization source is unique and not extended. It is 
applicable as long as the diffuse ionizing radiation field is not dominant and strongly inhomo¬ 
geneous. As examples of the capabilities of these new tools, we consider two very different 
theoretical cases. One is that of a high excitation planetary nebula that has an ellipsoidal shape 
with two polar density knots. The other one is that of a blister HII region, for which we have 
also constructed a spherical model (the spherical impostor) which has exactly the same H/1 
surface brightness distribution as the blister model and the same ionizing star. 

We present and comment line intensity maps corresponding to different viewing angles. 
We also use the computed line intensities to derive physical properties of the model in the 
same way as an observer would do for a real object. Eor example, we derive the “apparent” 
value of N/O for the entire nebulae and along spectral slits of different orientations. Eor this 
we take the electron temperature and density derived from the [N ii] 5755A/[N ii] 6583A and 
[O ii] 3726A/[0 ii] 3729A ratios, respectively, and we adopt the common recipe: N/0=N^/0^. 
Interestingly, we find that, in the case of our high excitation nebula, the derived N/O is within 
10-20% of the real value, even when the slit crosses the high density knots. On the other hand, 
for the blister HII region and its spherical impostor, we find that the apparent N/O is much 
smaller than the true one (about 0.68 and 0.5 of it, respectively). 

These two examples warn against preconceived ideas when interpreting spectroscopic 
and imaging data of HII regions and planetary nebulae. The tools NEBU_3D and VIS JS1EB3D, 
which will be made publicly available in the future, should facilitate the performance of nu¬ 
merical experiments, to yield a better understanding of the physics of aspherical ionized neb¬ 
ulae. 
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1 INTRODUCTION 


Most ionized nebulae have complex morphological structures. This 
is true not only for HII regions, which result from the ionization of 


molecular clouds b 

y newly 

born massive stars JKim & Koo 

|200i 

iReich etalJI 199(1 

O’Deli 

1200 it iMartm-Hemandez et alJ 

2003h. 

but also for planetary nebulae (see e.s. iManchado et al 

|l996|: 

iGomv et alJl999l:lMeixner et alj2004j). not to mention imaees ob- 


tained by the Hubble Space Telescope which show a large variety 
of complex structures. 

So far, almost all self-consistent physical analysis of ion¬ 
ized nebulae have used ID photoionization codes. This is due 
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to the fact that most photoionization codes are in ID. In many 
cases, with appropriate handling, these codes are actually able 
to infer important conclusions regarding the ionization structure, 
the heating or the elemental abundances of the ne bulae, even 
if the true geometrical structure i s more complex iQegg_eLalJ 
[l98j_|stasihska & SchaereijJjjQg GLuridiana & Peimberjl20(Ml 

iMoris^^^nTOO^Ijam^et all2004l). 

A few 3D photoionization codes now e xist jBaessgen_et_aL 

; iGruenwald et alJ Il997l : lErcolano et alJ ItOOsXTwoo^^I 
2004) . Quite naturally, the first 3D photoionization models applied 
to real nebulae have emphasized their success in reproducing se¬ 
lected line ratios, rather than their failure to reproduce all the ob¬ 
serva t ional c onstraints at the same time (e.g. iMonteiro et alJ200Cl 
l2004l2005l) . However, reproducing simultaneously all the obser- 
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vational constraints, and not only strong line ratios, is paramount 
to check our understanding of the physics of these objects. Fail¬ 
ure to reproduce let it be one line ratio, e.g. Fie ii/H/l or [O iii] 
4363/5007A may indicate that an important ingredient is missing 
e.g. in the description of the ionizing radiation field or in the en¬ 
ergy budget. Some recent work is starting to address issues related 
to the p hysics of ionized nebulae in which a 3D treatment is im- 

S it iErcolano et alJl2003cltt Iwood et al]l2004l : lErcolano et alJ 

There are however practical difficulties with 3D model fitting 
of real nebulae. One is the enormous parameter space to consider, 
especially in the case of complex morphologies. Another one is 
the fact that 3D photoionization models are very CPU-time and 
memory consuming: several hours are needed on clusters or super¬ 
computers to make a model (even if the time and memory can be 
reduced when modeling simpler geometries). Also, the difficulty in 
visualising the results does not help. 

It would therefore be useful to have at one’s disposal a quick 
3D photoionization code coupled with an efficient visualisation 
tool. This would allow one to explore more efficiently a complex 
parameter space, in the search of a satisfactory solution. This would 
also be of important tutorial value for experiments aimed at esti¬ 
mating the effects of the nebular geometry on the derived physical 
parameters (e. g. chemical composition, ionizing radiation field) or 
at educating oneself about how a 3D physical object may appear to 
an observer under various viewing angles. 

In this paper, we present a pseudo-3D photoionization code 
(NEBU_3D), which models aspherical planetary nebulae in a few 
minutes. The code uses various runs of a ID photoionization code 
and constructs from them a 3D photoionization model. The main 
limitation comes from the fact that the diffuse radiation is not 
treated e xactly, but in the ”outward-only” approximation as de¬ 
signed bv lTarten il967fi . where the diffuse radiation is assumed to 
be emitted outwards in a solid angle of In. Its associated visualisa¬ 
tion tool, VISNEB_3D, allows one to easily plot diagrams that are 
useful for understanding the physical conditions inside the models 
and diagrams that can be directly compared to observations. This 
tool can also be used to post-process the 3D grids produced by 
proper 3D photoionization codes. 

The code and the visualization tool are described in Section 2. 
The next sections are devoted to two applications: Section 3 treats 
the case of an ellipsoidal planetary nebula with two symmetrical 
polar knots and comments on the interpretation of observations of 
such knots in real nebulae in terms of chemical composition and 
excitation. Section 4 considers the case of a simple blister HII re¬ 
gion and compares it to a spherical HII region having the same Yip 
surface brightness distribution. Conclusions and perspectives are 
given in Section 5. 


2 NEBU_3D, A PSEUDO-3D PHOTOIONIZATION CODE 

NEBU_3D is a pseudo-3D photoionization code. It consists in 
a set of tools, developed under IDL (RSI), used to construct a 
3D model of an ionized nebula fro m a set of runs obtained by 
the ID photoionizatio n code NEBU jMorisset & Peauignolill99fit 
IPeauignot et alJEoOlll . It can easily be adapted to any other ID 
photoionization code. 

NEB U_3D is thus not a ful l 3D photoionization code like 
Mocassin Frcolano et alJ l2003all for example or the code by 
IWood et ^ J2004h . The main limitation of NEBU_3D is in the 
treatment of the diffuse field, which is not fully consistent. For ex¬ 


ample, NEBU_3D should not be used to model the ionization of 
shadows behind knots, as ID models cannot account in an accu¬ 
rate way for the ionization of the shadowed region by radiation 
produced by the surrounding material. On the other hand, when 
the radiation field is not changing drastically from one direction to 
the other, the use of a set of ID models does not lead to strong 
errors. The main adyantage of NEBU_3D is that it is faster by 3 
to 4 orders of magnitude than a full 3D code when running on 
the same computer. Note also that its predictions are more accu¬ 
rate th an those of the 3D photoio nization codes of iBaessgen et alJ 
h995h or iGruenwald et all il997h . The procedure to make a model 
with NEBU_3D is the following. We first define a cube in which 
we place the nebula. For each cell of the cube, we define the polar 
coordinates, according to the size of the nebula we are modeling. 

We run NEBU yarious times, changing the nebular density dis¬ 
tribution at each run according to an angular yariation law (NEBU 
also allows one to yary the chemical composition at any point in the 
nebula if one wishes to explore the effects of abundance inhomo¬ 
geneities). The set of outputs (such as the electron temperature, the 
emissiyities of selected lines at each radial position, etc...), which is 
defined according to the goal of the model, is then read and interpo¬ 
lated on the coordinate cube. In fact, if the nebula is axisymmetric, 
we only need to compute the interpolated yalues on l/8th of the 
cube, reconstructing the whole cube by rotating the results. With 
such a procedure, one can easily construct models of ellipsoidal or 
bipolar nebulae. One can also treat champagne flows or eyen irreg¬ 
ular nebular geometries proyided that the photoionization source is 
unique and that the effects of non-radial ionization are not expected 
to be important. 

The resulting emissivity cubes (one for each line of interest) 
can then be rotated and a projection on the sky plane is done to 
obtain surface brightness maps. 

We checke d that we can recoyer the whole set of the 108 mod¬ 
els published by lZhang & Kwol3 il998ll for the Yip images of ide¬ 
alized planetary nebulae, if adopting the same distribution for the 
gas. 


3 MODEL FOR AN ELLISPOIDAL PLANETARY 
NEBULA WITH TWO POLAR KNOTS 

3.1 Construction of the model and method of analysis 

A large fraction of planetary nebulae haye an ellipsoidal ap¬ 
pearance on the sky. Some of these, e.g. NGC 7009, NGC 
6826, NGC 3242, present symmetric knots or ansae. It has 
been proposed that in some cases the knot s may be nitrogen- 
rich tealick et ^JIl994l:lHyung & AlleJll 99.5h. although more re- 

cent analyses jAlexander & Balicldll997l : iGoncalyes et'^l2003t 

IPerinotto et all2004i) are much less categorical on this point. This 
requires deeper inyestigation, as the ionization correction factors 
may be different in the knots and in the surrounding material, due 
to different ionization conditions. The model presented in this sec¬ 
tion illustrates this issue. 

The density distribution is described by two components: a 
prolate ellispoidal nebula and polar knots. The ellipsoidal nebula 
is defined by an inner ellipsoidal cavity ri„{6) with an equatorial 
over polar axes ratio bja of 0.5. The density at the inner surface of 
the nebula is proportional to the inverse square of the distance to 
the ionizing star. Inside the nebula, the density decreases following 
a gaussian law. The polar knots are defined by adding a gaussian 
overdensity above the nebular background. 
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For a given polar angle 9 between 0 and njl {9 being 0 at the 
equator) the density distribution at the distance r from the ionizing 
source is then given hy: 

rm(9, r) = nl( [C^'^^\9, r) + r)) (1) 

\ nn(c') / 

where C'='*‘Pn6', r) and r) are dimensionless functions of the 

position describing the ellipsoidal and knot contrihutions to the 
density, defined as: 


C=“‘P^(e,r>r.„(0)) 


C='“P^(0,r<r,„(e)) 


r) 


exp - 


tsr{9) -3 
^ Strdh I 


0 


Crexp- 


Ar(9) - A-po,Ar 5 „(g) A 


exp - 


9-Till 

A'«_,.,e/(;r/2)j 


using: 

r-mid) = nn(0) - 

•^((hsin(6))2 + (acos(6))2) 
Ar{9) = r-r,M 


( 2 ) 

(3) 

(4) 


(5) 

(6) 


The relative intensity of the knot contribution is defined us¬ 
ing The size scales of the nebula and the knots are defined 
relatively to the thickness of an ionized shell Ar 5 ,r( 0 )* correspond¬ 
ing to a constant density nYi{9, ri„(9)) around an empty cavity of 
size r,„(6). Dimensionless coefficients are applied to these Arstr(9), 
namely Ssize which is used for the gaussian decrease of density in 
the ellipsoidal contribution, used for the position of the knot 
on the polar axes within the nebula and K^-size used for the radial 
size of the knot. Finally, A'a_size is the angular size of the knot, nor¬ 
malized to n/2. The reason for using such a density law is that a 
change in the physical parameters (nj^, ri„(0) and L.) whilst keep¬ 
ing the dimensionless coefficients constant allow us to change the 
ionization parameter keeping the global shape of the nebula. 

The model presented in this section is obtained using the 
following values for the physical parameters: = 2 10^ cm“^, 

An(0) = 10'’ cm, and the following values for the dimensionless 
parameters: a = 2, b = 1, = 5.0, S^ize = 1-7, A'pos = 0.3, 

Kz-size = 0.1, and Kg-^ize - 0.078. 

We have adopted the following chemical composition for the 
entire nebula, including the knots: Ft:Fte:C:N:0:Ne:S:Ar= 1:0.1:3 
10-*:1 10-‘*:6 IQ-^LS 10-‘':1.5 10-^:3 10-^ by number. 

The characteristics of the ionizing star are T^f f = 80 kK and 
L, = 10’® erg/s. 

We have run 40 models with NEBU to cover the domain 
[0, ;r/2] for 9. The number of steps for the NEBU runs is between 
60 and 85, depending on the density distribution associated with 
each angle. 

The results of the NEBU runs are then interpolated in the co¬ 
ordinate cube. The cube in which the interpolation id done is com¬ 
posed of 200’ cells. For the nebula considered here, this correspond 
to a size for each cell of 2.28 10'® cm. The variables that are in¬ 
terpolated are: the electron temperature and density, the hydrogen 


In the case of a constant density shell around a cavity, we have: Arstr = 


finii- 


36(1 


■ + 1)'” - 1) 
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Figure 1. Cut in the hydrogen density distribution for the model presented 
in Sect. 1.3. li The solid line indicates the recombination front of the nebula. 
The grey scale is in cm“’. 


density, the ionic abundances of N'^ and O^, and the emissivities of 
the line of interest, namely Fly0, He ii 4686A, [N ii] 5755A, [N ii] 
6583A, [O It] 3726A, [O ii] 3729A and [O iii] 5007A. 

Figure Q is a 2D map of the hydrogen density distribution of 
our model in a plane containing the polar axis. The star is at coor¬ 
dinates (0,0), and the recombination front is also plotted. Figure|5| 
is a 3D isodensity representation of the model. The chosen density 
level allows one to clearly distinguish the main body of the nebula 
and the polar knots. 

After having computed the model, we generate line intensity 
maps corresponding to “pole-on” and “face-on” viewing angles, by 
summing up the corresponding emissions on each line of sight. An 
example of such “pole-on” and “face-on” contour maps is shown 
in Fig.|^ Line ratio maps are then obtained by performing a pixel 
by pixel division of the two line intensity maps. 

We then consider the model as if it were a real nebula, and 
apply to it the standard abundance analysis using the same atomic 
data as in NEBU. Eor each pixel of the [N ii] 5755/6583A map we 
compute the “apparent” electron temperature and for each pixel of 
the [O It] 3726/29A we compute the “apparent” electron density. 
These values are used to compute for each line of sight the “ap¬ 
parent” N^/0^ ratio from the [N ii] 6583A/ [O ii] 3726,29A maps, 
which is identified to N/O, following the common practice. 


3.2 The results 

3.2.1 Synthetic nebular images 

Fig .|^shows intensity maps of the model nebula seen in a direction 
perpendicular to the polar axis. Four emission lines are represented: 
Bj3, He II 4686A, [O ii] 3726,29A, and [O in] 5007A. Such a pre¬ 
sentation allows one to distinguish immediately various character¬ 
istics of the model images. For example, it is clearly seen that the 
H/3 and [O ii] 3726,29A images are more extended than the He ii 
4686A one, and that the extension of the [O iii] 5007A is interme¬ 
diate. One also sees that the emission is stronger along the equator. 
This is due to the density being smaller at high polar angle. The in¬ 
ner “egg-shape” seen in H/3 and [O iii] 5007A results from the elip- 
soidal inner empty cavity. The polar knots are well visible in the H/3 
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Figure 2. 3D representation of the nebula. An isodensity surface is drawn, 
showing the equatorial density enhancement and the two polar knots. On 
the faces we have represented the H/3 surface brigthness contours for 2 ori¬ 
entations of the nebula: on the left side for the nebula seen with the polar 
axis paralle to the sky, on the top for the nebula seen pole-on. The slits used 
to determine pai'ameters from emission line ratios are indicated. 


map and even more conspicuous in the [O ii] 3726,29A map (the 
same applies for the [N ii] 6583A line map, not shown here), due to 
enhanced recombination in the knots. Even if the model presented 
in this section does not attempt to reproduce a “real” nebula, such 
emission enhancement of the low ionization speci es can be seen in 
real n ebulae showing ANSAE such as NGC 7009 iOoncalves et alJ 
l2003h . 

3.2.2 Local properties 

Figs. El and 13 show the variation of some local properties in the 
model grid {T^, n„ (N+/0+)/(N/0), N+/N, and 0+/0) along the 
polar and equatorial axes respectively. We see that the maximum 
density of the knot is larger by a factor about 3 than the inner den¬ 
sity along the polar axis (Fig.|4j and only sligthly larger than the 
inner density on the equatorial plane (compare the dotted line in 
Figs.|3and|3. Note that the zone of maximum density in the knots 
lies out of the ionized region. The electron temperature does not 
show significant structure apart from the natural drop in the recom¬ 
bination zone (see the solid line in Figs.|3and|3. The ionization 
fractions N^/N (dashed-dot line) and 0^/0 (dot-dot-dashed line) 
increase away from the central star, as naturally expected, with 
the increase being more important in the zone of the dense knot. 
The (N^/0^)/(N/0) ratio (dashed line), which is generally assumed 
equal to one in abundance analysis, is actually never strictly equal 
to one in our model. It is of the order of 0.6 - 0.7 in the high ion¬ 
ization zone and increases to values of about 2.5 on the back side 



Figure 3. Synthetic intensity maps of the ellipsoidal nebula when seen per¬ 
pendicular to the polar axis. Each quadrant corresponds to an emission line: 
Upper left: Hfl, upper right: He ii 4686A, lower left: [O ii] 3726,29A, and 
lower right: [0 iii] 5007A. The intensity maps are normalized so that the 
maximum for each line is the same. 



pixels 


Figure 4. Variation of some local properties in the 3D grid of the model, 
along the large axis of the nebula (i.e. crossing the knots). Solid line: 
(10‘*K), dotted: «h (10^ cm“^), dashed: (N+/0^)/(N/0), dashed-dot: N+/N, 
and dot-dot-dashed: 0*10. Note that in this figure the pixels range from 80 
to 120. 


of the knot as well as in the ionization front. In the latter case, the 
electron temperature is already low, so that this increase will have 
no effect on observed line ratios. The strong maximum at the back 
side of the knot is probably exagerated in our model, since the dif¬ 
fuse ionizing radiation is not treated properly. 
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pixels 

Figure 5. Same as Fig. |3]but along a small axis of the nebula. Note that in 
this figure the pixels range from 40 to 90. 


3.2.3 Line of sight intensities 

Figs.|^to|^show how the above local properties translate into the 
intensities of emission lines as would be observed with a narrow 
slit. Fig.l^corresponds to the case where the slit is on the large axis 
and crosses the knots (slit SI in Fig.|3- Fig.Qcorresponds to the 
case where the slit is along the small axis (slit S2 in Fig.|2j. Fig.|^ 
corresponds to the case where the nebula is seen pole-on with a slit 
passing through the center of the knots (slit S3 in Fig.|5J. The in¬ 
tensity axes are in arbitrary units, but the relative line intensities are 
meaningful (applying the coefficients given in the caption). Fig.|^ 
shows that all the line of sight intensities present an enhancement 
at the position of the knot. However, the enhancement is moderate 
in the case of [O iii] 5007A (dot-dashed line) due to the opposite 
effects of density enhancement and recombination of in the 
knot. In the case of H/1 (solid line), the intensity is roughly propor¬ 
tional to the integral of the squared density along the line of sight, 
so the enhancement is larger than for [O iii] 5007A (with respect 
to the nearby lines of sight, the enhancement is of about a factor 3 
while it is only of 50% for [O iii] 5007A). The intensities of [N ii] 
6583A (dashed line) and [O ii] 3726,29A (dotted line) are more 
strongly enhanced (by about a factor 4-5) because the N+/N and 
0^/0 ratios are larger in the knot than in the surrounding mate¬ 
rial. The inner cavity and the equatorial density enhancement are 
clearly seen in spectra along the small axis and when the nebula is 
seen pole-on (see Figs.0and|^. On the other hand, they are not 
seen when the slit is along the large axis (Fig.|^, because the lack 
of emission in the cavity is compensated by the stronger emissivity 
in the equatorial plane where the gas is overall denser than on the 
polar axis. Note that in our model, when the nebula is seen pole-on 
(Fig-ID, the bump corresponding to the knot is barely seen and the 
emission is dominated by photons arising from dense regions from 
the body of the nebula for which the emission measure is large in 
such a configuration. 


3.2.4 Apparent properties 

Fiss.l^to ll ll show the values of (solid line), Ne (dotted line) and 
N/O (dashed line) determined using the line of sight intensities of 
[N It] 5155k, [N It] 6583A, [O ii] 3726A and [O ii] 3129k lines in 
Figs. |^to|^ respectively. The procedure to derive these quantities 



pixels 

Figure 6. Variation of the intensity of some lines along a nan'ow slit ori¬ 
ented from the ionizing star to the polar knot (slit SI in Fig.l^. Solid: lOx 
H/?, dotted: 3x 10 ii] 3726,29A, dashed: 3x [N ii] 6583A and dot-dashed: 
10 III] 5007A. The knot is in the vicinity of pixel 110 
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Figure 7. Same as Fig.|^ but through a narrow slit in the equatorial plane 
(slit S2 in Fig.|2j. 

has been described in Sect. B.ll and corresponds to what an observer 
would do using the line intensities from long slit spectra of a real 
object. We call these quantities “apparent” properties of the model. 

As seen in Fig.|3 in the case of a slit along the polar axis, the 
apparent electron density decreases steadily from the center until 
the position of the knot is reached where a sharp peak is seen. On 
the other hand, when the nebula is observed through a slit aligned 
on the small axis (Fig. GHJ. the apparent density is constant over 
most of the nebular apparent surface except the outermost zone 
where it decreases abruptly. The presence of the central cavity, of 
course, cannot be guessed from such line ratios, and can only be 
inferred from the radial distribution of line fluxes such as presented 
in Figs.|^to[^ The fact that the apparent density is uniform if the 
nebula is seen along a small axis is simply due to the fact that most 
of the [O ii] 3726,29A emission comes from a thin external zone. 
On the other hand, along the polar axis, the apparent density in the 
central parts is strongly affected by the dense equatorial zone, and 
therefore the apparent density decreases outwards. 
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Figure 8. Same as Fig.|^ but for the nebula seen “pole-on” through a naiTow Figure 10. Same as Fig.|^ but for a a narrow slit oriented along the small 

slit passing through the center of the knots, (slit S3 in Fig.|3. The knot is axis (slit S2 in Fig.l^. 

centered at pixel 0. 



Figure 9. Apparent variation of some properties derived from the emission 
line ratios along a narrow slit oriented along the polar axis and crossing the 
knot (slit SI in Fig.|3. Solid line: (10“*K), dotted: (10^ cm“^), dashed: 
N/0 normalized to the “true” value of 1/6. 


The apparent N/0 value (obtained assuming N/0=N+/0^) is 
never identical to the true value. As seen in Figs. and 1101 the 
ratio (apparent N/O) / (true N/O), i.e. the dashed lines is 0.85 in the 
whole nebula, except close to the knot, where it drops to 0.78 before 
increasing up to a value of 1.2 at the far side of the knot. This means 
that from the observations of the emission lines shown in Fig.|^ one 
would infer a modest variation (up to 40% in this model) of the N/O 
ratio across the nebula. Taking into account the slit width and the 
seeing, an observer would probably not distinguish any abundance 
variation in such a case. Notice that the apparent increase of N/O at 
the edge of the nebula to values up to 2.5 times the true value (pixel 
numbers greater than 160 and 80 in Figs.Mand llOl resp.) are not 
relevant, as the emission of all the lines is very low at this position 
in the nebula (see Figs.|^andQresp.), and such an increase could 
not be seen. 

Fig . ll2l shows the map of the bias in N/O for the case of a neb¬ 
ula whose polar axis is in the plane of the sky. The pixels where the 
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Figure 11. Same as Fig.|^ but for a “pole-on” nebula (slit S3 in Fig.|3. 


intensity of [O ii] 3726,29A is lower than 0.05 times its maximum 
value are set to black, to avoid misinterpretation of the increase of 
N/O at the edge of the nebula, where no emission arises. Similarly 
to what is seen in Fig.^| we see the apparent increase of N/O just 
behind the knot. 

Fig. ll3l 2D histogramm of the number of cells with a given 
value of Ue n(N^) and n(O^), for all the cells in the nebula where 
0^/0 > 0.02. The vertical scale on the right gives the correspon- 
dance between the grey level and the log of the number of cells. The 
parts of the nebula emitting strongly in [O ii] 3726,29A and [N ii] 
6583A are far from the (0,0) origin, the parts emitting weakly in 
these lines are close to it. If the assumption that N/0=N+/0^ were 
true, then all the points would lie on the y=x axis. That is clearly 
not the case, the majority of points are either under or above the 
first bisector. 

Note that the N/O derived for the entire nebula from an in¬ 
tegrated spectrum is 0.93 in this case, which happens to be very 
close to the true value, due to counterbalancing effects in the vari¬ 
ous zones. This example shows that, even if N^ and are present 
in small proportion in a high excitation nebula, the overall N/O 
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Figure 12. Image of (apparent N/0) / (true N/0) obtained from the emission 
line maps of Fig. 1. 
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Figure 13. 2D histogramm of the number of cells with a given value of 
(N^/N) He and (0^/0) n^, for all the cells in the nebula where 0^/0 > 0.02. 
The vertical scale on the right gives the con'espondance between the grey 
level and the log of the number of cells. 


abundance ratio derived by classical empirical methods may be rea¬ 
sonably accurate. 

Whether our results on the apparent N/0 value are valid in 
general (or at least for the majority of planetary nebulae) would 
require a scan of the entire parameter space. 


4 MODELS FOR A BLISTER HII REGION AND ITS 
SPHERICAL IMPOSTOR 

The idea that HII regions may be described by a champagne flow 
model has been widespread since the paper hv IZnckermanI il97.3ll 
who showed that the Orion nebula is not a sphe re but has the 
structu re of a blister, and the theoretical paper bv iTenorio-Taglel 
<1979h , who showed how such a structure is likely to arise. How¬ 
ever, very little photoionization work for such geome tries has been 
perfor med. Among the few exceptions aretheworkof| Binetteet all 
( I 2 OO 3 I) on the Orion nebula, and that of iMorisset et alJ i2002h on 
G29.96-0.02 to date. 


As a matter of fact, most modelling is still done either in plane 
parallel approximation, or in spherical symmetry. Here, we con¬ 
struct a photoionization model for a blister HII region ionized by 
a single star. The Hy6 image of such a model when seen face on is 
circular, and the object could be mistaken for a sphere. 

In this section, we discuss the properties of a blister model 
and compare them with those of a spherical model that has the same 
ionizing star, the same chemical composition and the same total Hj8 
luminosity and surface brightness distribution. We call this second 
model the spherical impostor. 

4.1 The models 

The blister model presented here is a plane parallel slab, with 
density increasing e xponentially wi t h the distance to the ioniz¬ 
ing star, as found by IWen & O’Delll il99.5ll for the Orion nebula: 
”h(z) = "h exp((z - ZbilQ, where z is the normal distance to 
the star and Zi,T are parameters taken to be equal to 2 10'* and 5 
lO'^ cm respectively. The value of nj^ in our model is 3 10'* cm“^. 
The star is assumed to radiate as a blackbody with an effective tem¬ 
perature of Te// = 37000 K and a luminosity of 3.5 10^*erg/s. 
The chemical composition is solar. The gas density distribution 
law of the spherical impostor has been obtained by trial and er¬ 
ror and is the following: n^(r) = exp(-(r/C)^) x B~l{d + D^), 

where = 2.610* cm“*, B = 1.5 10** cm, C=1.2 lO'* cm, andD=6 
10 '* cm. An empty inner cavity around the ionizing star of radius 
10'^ cm is assumed. 


4.1.1 Spatial distribution of the line intensities 

Figure lT^l presents line intensity maps of the models, in a way simi¬ 
lar to Fig.0 but with He ii 4686A replaced by He 1 6678A, since the 
stellar temperature is only of 37000 K. The left panel corresponds 
to the blister model, seen tangentially to the slab. In this panel the 
[O ii] 3726,29A and [O in] 5007A maps are upside down, for a 
better comparison. Note that the blister model is not a plane par¬ 
allel model, since the radiation comes from a star whose distance 
to the slab is not infinite. We see that, for each line, the surface 
brightness increases towards the ionization front. This is due to the 
adopted density law, which increases with distance to the star. We 
also see that the surface brightness close to the ionization front is 
larger in the direction perpendicular to the slab, and decreases with 
increasing polar angle. This is due to the dilution of the radiation 
which is more important at higher polar angles, so that the ioniza¬ 
tion front occurs in a less dense zone. As expected, the extensions 
of the [O III] 5007A and He 1 6678A emitting zones are smaller than 
those of the [O ii] 3726,29A and emitting zones. This is better 
seen in the maps corresponding to the face-on blister, which are 
presented in the central panel of Fig. ESI The right panel of Fig. ESI 
presents the line intensity maps of the spherical impostor model. By 
construction, the map is identical to that of the face-on blister. 

The maps for the remaining lines are significantly different. 
They are better discussed by considering also Fig.^] which com¬ 
pares the HJ3, [O ii] 3726,29A and [O iii] 5007A line intensity vari¬ 
ations for the face-on blister (dotted lines) and for the spherical im¬ 
postor (solid line). We seen that the excitation level of the spherical 
impostor is much higher than that of the blister. The reason is that, 
in the spherical impostor, the dense matter is close to the star and 
receives a strong radiation field, while in the blister, the dense mat¬ 
ter is far from the star. In the spherical impostor, the [O iii] 5007A 
intensity decreases rapidly away from the central line of sight, be¬ 
cause the ionizing photons are quickly exhausted due to the 
























8 C. Morisset et al. 



Figure 14. Synthetic H/3, He i 6678A, [O ii] 3726,29A and [0 in] 5007A surface brightness maps from upper left to lower right. Left panel: the blister model 
seen tangentially (in this case the maps for [O ii] 3726,29A and [O iii] 5007A are upside down); central panel: the blister model seen face on; right panel: the 
spherical impostor model. 


dense inner region, while in the face-on blister, the variation of the 
ionizing radiation field penetrating the gas is much smoother. The 
[O ii] 3726,29A map also shows conspicuous dilferences between 
the two cases^. In the spherical impostor, the curves related to the 
[O ii] 3726,29A intensity first gently increases away from the cen¬ 
tral line of sight, due to the radial increase of the population 
convolved with the decrease of the gas emission measure on the 
line of sight (the integral of the square of the gas density). Farther 
away, all the oxygen is in the form, and the [O ii] 3726,29A 
intensity decreases. On the other hand, for the face-on blister, the 
[O ii] 3726,29A intensity decreases steadily outwards. Indeed, the 
effect of the emission measure decreasing with increasing polar an¬ 
gle dominates over the ionization effect. This result is of course 
strongly dependent on the gas density law taken for the slab. With 
a much weaker density increase, one might observe a behaviour 
qualitatively similar to the one shown by the solid line in the right 
panel of Fig. 11.51 

Figure fT^ is similar to Fig.^| but shows line ratios instead of 
line intensities: from left to right, He 1 5876A/II;6, [O iii] 5OO7A/II;0 
and [O ii] 3726,29A/Flj8. We see that [O iii] 5OO7A/II;0 decreases 
abruptly for the spherical impostor, while it decreases very slowly 
for the blister. Similarly, [O ii] 3726,29A/II/3 increases abruptly for 
the spherical impostor, and slowly for the blister. Indeed, the ion¬ 
ization parameters of the sampled zones are much higher in the 
spherical impostor than in the face-on blister, resulting in a much 
higher central [O iii] 5007A/F1/1 ratio and a sharper ionization strat¬ 
ification. As for the He i 5876A/Hj8 ratio, it decreases outwards, as 
expected. In the spherical impostor, He i 5876A/H;6 abruptly goes 
down to zero at about 90 pixels, while the decrease is much slower 
in the face-on blister, again because of the sharper ionization strati¬ 
fication in the spherical impostor. The total He i 5876A/H;S value is 
0.029 for the spherical impostor , i.e. slightly larger than the value 
of 0.027 for the blister. 

This example shows that, if one were to construct a photoion¬ 
ization model for a blister HII region seen face on, and assumed a 
spherical geometry, one would run into severe difficulties matching 
the observed line ratios as a function of projected distance to the 
star. 


^ In the case of the blister model, the [O ii] 3726,29A lines are not smooth 
for numerical reasons, the size of the region emitting these lines beeing of 
the order of a few cells (See Fig. EH. We checked that the values are close 
to the outputs of the original NEBU runs. 



(OVO).n, 


Figure 18. Same as Fig. [HI but for the blister model. 


4.1.2 Comparison of apparent properties 

We now compare the apparent physical properties of the models 
(Te, We, and N/0) derived from empirical methods in a way simi¬ 
lar to what we presented in Sect. l3.ll Figure fTTl shows the electron 
density derived from [O ii] 3726/29A (left panel), the electron tem¬ 
perature derived from [N ii] 5755/6583A (middle panel) and the 
empirically derived N/O ratio (right panel). In both models, we ob¬ 
serve a general decrease of the apparent electron density away from 
the star. In the case of the blister, this apparent density decrease is 
due to the fact that, as the polar angle increases, the regions that 
are sampled by the [O ii] 3726,29A lines are of decreasing density, 
due to the increasing dilution of the ionizing radiation field. Note 
that the apparent densit y decrease we find fo r our model is very 
similar to that found bv lWen & O’Delll il99.5h for the Orion neb¬ 
ula. In the case of the spherical impostor, the true hydrogen density 
decreases away from the central star, while the apparent density de¬ 
duced from the [O ii] 3726,29A line ratio is constant within the first 
70 pixels, because the emission lines are strongly weighted by the 
inner spherical skin of the zone. 

The observed apparent electron temperature is remarkably 
similar in both models within the first 100 pixels. This however 
may be fortuitous, and might not occur for different density distri¬ 
butions, because the electron temperature depends on the relative 
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Figure 15. Emission line intensity variations (from left ro right: Hj3, [O iii] 5007A, [O ii] 3726,29A). Solid line: spherical model, dashed line: blister model. 





Figure 16. Emission line intensity valuations, normalized to (from left ro right: He i 5876A, [O iii] 5007A, [O ii] 3726,29A). Solid line: spherical model, 
dashed line: blister model. 
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Figure 19. Same as Fig. ini but for the spherical impostor model. 


proportions of different ions on the line of sight. In the outer parts, 
on the other hand, there is a conspicuous difference. While the ap¬ 
parent temperature continues to slowly decrease in the case of the 
blister, it significantly increases in the spherical impostor. This is 
due to a more important hardening of the ionizing radiation field in 
the latter case. 

Finally, the derived N/O ratio is equal to the input value only in 
the outer part of the spherical impostor, where the flux is very small. 
In the brightest regions of both models, the apparent N/O ratio is 


smaller than the true one by significant amounts (about 0.25 and 
0.65 times the true value in the spherical and blister models, respec¬ 
tively). The apparent N/O ratio for the entire nebulae is 0.68 (0.50) 
times smaller than the true one in the case of the blister (spherical 
impostor). 

Figures fT^ and fT^ show the 2D histogram of the number of 
cells with a given value of «e n(N+) and n(0*] within the neb¬ 
ula (blister and spherical impostor resp.), in the same way as in 
Fig. 11.31 In both cases, an important amount of the gas is in a sit¬ 
uation where N+/0^=N/0 (points on the first bisector), but a non 
negligible amount is also located far from this line, causing the dis¬ 
crepancy between the “apparent” and true values for N/O. Note that 
in the case of the impostor, the fact that the model is spherical leads 
to a very narrow line in the 2D histogram. 

The main conclusion of this section is that the apparent proper¬ 
ties of two models that differ only by their geometries and in such 
a way to conserve the lift surface brightness distribution may be 
significantly different. We also note that, in the cases we consid¬ 
ered, the N/O ratio derived for the entire nebulae using the clas¬ 
sical recipe that N/0=N+/0^ is significantly below the true value, 
even though N^/N and 0^/0 are not very small. In the present case, 
this is due to the spectral distribution of the ionizing radiation field 
(with realistic stellar atmospheres, the derived N/O would likely be 
different). 


5 CONCLUSION AND FUTURE WORK 

We have presented a quick pseudo-3D photoionization code, 
NEBU_3D, and its associated visualization tool, VIS_NEB3D. This 
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Figure 17. Vaiiation of some “apparent” parameters (from left ro right: n^, T^, and N/0), determined from the emission line ratios shown in Fig. 1151 Solid 
line: spherical model, dashed line: blister model. 


code is able to treat a large variety of problems in a time that is 
by orders of magnitude faster than the time required by proper 3D 
photoionization codes which treat the diffuse radiation in an exact 
manner. NEBU_3D is not able to correctly treat cases where the dif¬ 
fuse radiation field is important, such as shadows behind optically 
think clumps or very complex geometries. On the other hand, it is 
more accurate than 3D photoionization codes that work under the 
on-the-spot approximation for the diffuse ionizing radiation. 

The main virtue of NEBU_3D over accurate 3D photoioniza¬ 
tion codes is that it allows an efficient exploration of a large pa¬ 
rameter space implied by non spherical geometries. This should be 
useful as a first approach for a tailored 3D photoionization mod¬ 
elling of real nebulae. Also, it allows one to easily run models 
with various geometries, and to get educated on the properties of 
such models which are not easy to predict without computations 
and therefore to exercise more care when interpreting observational 
data. 

We have treated two examples to illustrate the possible ap¬ 
plications of our code. One is the case of an ellipsoidal planetary 
nebula with two knots symmetrically located on the polar axis. We 
have shown and commented line intensity maps corresponding to 
different viewing angles. We have also used the computed line in¬ 
tensities to derive physical properties of the model in the same way 
as an observer would do for a real object. For example, we have 
derived an apparent value of N/O both for the entire nebula and 
along spectral slits of different orientations, assuming the common 
recipe that N/0=N^/0^. Interestingly, we found that, for this high 
excitation planetary nebula, for which one would expect the N/O to 
be less accurate than for lower excitation objects, the apparent N/O 
value is very close to the true one for all the observational set-ups. 

The other example is that of a blister HII region, and of its 
spherical impostor, i.e. a spherical model which has the same H/l 
flux surface brightness distribution, and the same ionizing star. We 
have shown that the properties of these two models differ drastically 
in many respects. We have also computed the apparent N/O ratios 
and found that they are significantly smaller (by factors 0.68 and 
0.50 for the blister and the spherical impostor resp.) than the true 
values. 

These two examples should serve as a warning against pre¬ 
conceived ideas when interpreting spectroscopic and imaging data 
of planetary nebulae and HII regions. 

The velocity fields of a blister HII region and of a spherical 
HII regions are however expected to be quite different. They should 
allow one to discriminate various models. We will address this issue 
in a future paper. 


We plan to make the NEBU_3D and VISNEB_3D codes pub¬ 
licly available in due time. 
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